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ABSTRACT 


Previous quantitative results obtained at M.I.T. on the 
bending of steel by laser line heating to produce simple 
angular deflections have been verified and expanded. 
Deflections and residual stresses in compound shapes formed 
by laser line heating of steel are measured and analyzed. 
From the meaSured deflections, algorithms are developed 
allowing the construction of idealized compound shapes 
including the dish, saddle, cone, and sine. All algorithms 


are developed from Single pass laser line heating data. 
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CHAPTER ONE 


INTRODUCTION 


It has been substantially demonstrated by Masubuchi [1], 
McCarthy [2], et. al. that bending of steel plates by laser 
line heating is practicable for simple one dimension bending 
and possible for bending complex shapes. In some ways laser 
bending is preferable to mechanical forming. These include: 
(1) Minimal material degradation [3],[4]; 

(2) The ability to determine deflection angles with heat in- 
put parameters; 

(3) The ability to automate the process fully with a closed- 
loop feedback control system. 

References [1] and [2] also show that complex shapes such as 
the dish, saddle, Sine, and cone which are not readily formed 
as a unit by mechanical processing, can easily be formed by 
laser line heating. Much literature exists extolling the 
advantages of line heating over mechanical forming for 
complex shapes. However, to date most attention has been 
focussed on the use of flames for bending steels. 

Historically, line heating by flame has been an art 
rather than a science. Attempts have been made to quantify 
the physics of flame bending including the description of 
plastic and elastic phenomena by boundary integral equations 


and finite element methods [5]. In spite of the application 





of advanced and sophisticated mathematical tools to the field 
of flame and laser bending, Simple and easily programmable 
algorithms to predict the curvature of complex shapes by 
flame or laser bending do not exist. 

Lasers have been used for decades to produce an intense 
and concentrated line or point heat for welding or producing 
micro and mini-holes. The technology already exists to con- 
trol accurately the heat input and heating pattern of a laser 
by optically defocussing the beam. This high degree of 
control indicates the major advantage of laser line heating 
over flame line heating. Specifically, a laser system is, by 
its very nature, quite adaptable to a closed-loop control 
System for thermal bending of steel into complex shapes 
including those with well-defined inflection points. 

Today 1s a period of tranSition in manufacturing from 
Simple open-loop and closed-loop control systems to the 
Gevelopment and implementation of "expert" and robotic 
Systems. The ability to control and to predict accurately 
relevant thermal parameters is critical if plate bending by 
line heating iS to become an accepted automated manufacturing 
process, 

In Spite of the expense of a high power laser system, its 
use in industry 1S increasing rapidly. Discussions with one 
major high power CO, industrial scale laser manufacturer 
revealed that, for that company, the sale of large scale 


industrial lasers iS increaSing about thirty per cent per 
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year. Many companies who would not be willing or able to 
apply the large capital expenditures required for a laser 
system might invest in an automated flame bending system. 
For that reason, the algorithms developed in thiS paper were 
designed to be applicable to flame bending with obvious 
adjustment to power-deflection parameters. 

The purpose of the research discussed in thiS paper is 
threefold. First, to prove the reproducability of laser line 
heating by predicting and then reproducing one of the experi- 
ments discussed in references [1] and [2]. Secondly, to 
measure the residual stresses found in one of the Shapes 
previously formed by laser line heating. Finally, and most 
importantly, to develop simple, linear algorithms which 
utilyzes the extenSive, Simple, one dimensional data base to 
mathematically construct complex shapes. Algorithms for the 
dish, saddle, cone, and Sine shape are separately developed. 
They are empirically formulated from shapes previously formed 
by laser line heating. These algorithms are then used to 
predict the deflections of shapes to be formed in subsequent 


experiments. 
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CHAPTER TWO 


PREVIOUS AND FUTURE RESEARCH 


2.1 Long-Term Research Goals 

A discussion of future research is normally included at 
the end of a paper. In this case, in order to understand 
adequately the intent and scope of current research in laser 
forming of steel plates, one must be aware of the long-term 
research goals of which the work reflected in thiS paper 1S 
but one step in the process. 

The ultimate goal of current research in laser bending of 
steel plates is the development of a closed-loop and/or 
expert system to fully automate the process of forming simple 
and complex shapes frequently found in the shipbuilding 
industry. Although the focus of research has been on uSing a 
high power (6-15 kW) laser as the line heater, it iS antici- 
pated that the method will be readily adaptable to flame 
bending. Several subsystems must be developed in parallel in 
Order to accomplish this: 

a) An extenSive data base for one dimensional bending must be 
established correlating plate deflections to heating and 
material parameters. 


b) Algorithms must be developed so that the system has 


initial values and heating patterns from which to start and 
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from which the system could estimate correct curvatures for 
complex shapes when encountering new sets of requirements. 

c) A "real-time" deflection measurement subsystem must be 
Geveloped to give closure to the feedback system. 

Gd) Computer programs and methodologies must be developed to 
control and interface the database, algorithms, and the real- 
time deflection meaSurements. Also, they must be capable of 
interfacing with the bending laser and movable bedplate 
systems. 

e) Prototype total system development and assembly. 
References [1] and [2] and part of this paper show that step 
"a" has been accomplished. A low power laser/video camera 
interferometer system 1s being developed at this time. It 
Should accomplish step "c",. The main purpose of this paper 
1s to describe the effort to accomplish step "b"; namely, the 
Gevelopment of the algorithms. Subsystem "d" and total 
System integration will follow from past and current 


research. 


2.2 Review of Previous Research 
References [1] and [2] discuss in detail the results of 
experiments on the bending of steel plates by a high power 


CO5 laser. Figure 2.2-1 is a sketch of the basic method of 


laser bending a plate. 
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Figure 2.2-1: Diagram of a Basic Laser Bending Apparatus 
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The steel plate restS On a movable bed plate such that the 
stationary beam irradiates the surface of the steel as a line 
pass. Primary deflection 1S perpendicular to the laser line 
Basse the resulting deflection angle, @-, is the sum of 6] 
and @5,. All deflection angles given in [1] and [2] are in 
terms of @¢- 

Many factors can affect the amount of curvature. These 
include: 
1) Plate thickness and type of steel. 
2) Heat input which is, itself, dependent on the laser power, 
optical path energy loss, absorbance fraction by the plate of 
the delivered infrared power, and the speed of the plate 
underneath the laser beam. 
3) Plate restraints such as clamping and free edge effects. 
4) Plate cooldown rate. 
5) Laser beam spot geometry may have an effect by changing 
the power density delivered to the plates surface. 
Since metals are highly reflective to the infrared radiation 
(~10 microns) emitted from a CO laser, the plate surface 
must be coated with an absorbing substance. Black spray 
paint is currently being used for that purpose. Optical 
losses are approximately ten per cent for the 15 kW laser 
used at the Naval Research Laboratory. 

Heat input is being represented by the following 


parameters: 
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P/A/V and P/twW 

where P is the laser power is kWatts less optical loss, V is 
the plate speed in inches-per-minute andt is the thickness 
of the plate in inches. The deflection angle, O¢, 1S 
correlated to specified P/W and/or P/tW. The deflection 
angle has been measured by dial gages and a laser 
interferometer ([1],[2]). Table 2.2-1 and figures 2.2-2 and 
2.2-3 are extracted from reference [1]. Clearly, the results 
represented in the graphs and table can be programmed into a 
computer with routines for estimating deflections at 
parameters not previously found by experiment. These results 
with the additional data measured and described later in this 
paper form the data base from which the algorithms for 
complex shapes are developed. 

The free edge of the plate where residual stresses must 
vanish have been shown ([1],{2]) to have a pronounced effect 
on the magnitude of the deflection angle for a given set of 
parameters. Distortion angle decreases relative to that 
obtained on the interior of the plate at the edges of a 
plate. However, approximately one beam diameter (1.5 inches, 
nominal) away from the edge, this effect is small and will 
generally be neglected in this paper. For purposes of the 
analyses in this paper the effects of the cooldown rate, edge 
restraint, and beam spot pattern also will be neglected 


unless otherwise specified. 
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Table 2.2-l: Process Parameters and Test Results Obtained 
in Parametric Studies 
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CHAPTER THREE 


RESIDUAL STRESSES IN COMPOUND SHAPES 
FORMED BY LASER LINE HEATING 


3.1 Introduction 

It 1S a well-known phenomenon that any forming of a 
metal will impose residual stress in that metal. The failure 
of a component is often abetted if not caused directly bv 
residual stresses. When no external load is applied nor any 
external restraint 1S present residual stresses can still be 
present. One cause of these residual stresses is a 
differential temperature in the metal. It 1s believed that 
the deflections resulting from line heating are caused by a 
combination of plastic and elastic behavior. Since residual 
Stresses are limited by the yield strength of the material, a 
measurement of the residual stresses in a compound shape 
Should indicate the relative contribution of elastic strain 
to the total deflection. Furthermore, for critical 
components where failure would cause the loss of the system, 
minimizing or removing residual stresses might be required. 
This 1S particularly true if the compound shapes are to be 
welded. 
3.2 Theoretical Relationships 

Metals are susceptible to residual stresses caused by 
uneven plastic strains. Masubuchi [6] developed some 
fundamental relationships for analyzing residual stresses in 


welds. Assuming perpendicular stresses are not significant 
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the following development should apply to line heating. he 
should be recognized, however, that since the heating from a 
laser is localized to the surface of the plate, some through 
thickness residual stresses probably exist. The following 
derivation 1s taken from reference [6]. 

a. Total strain resultS from elastic and plastic 


components: 
Ex = Ex(e) + Ex (p) 
Sy Se) vermy,(p.) 


yxy = Yxyle) + Xky(p) 


where €y, Ey rxy are components of the total Strain 


Ex(e)*s Ey(e)s Xxy(e) are components of the elas- 
mire sit Gain 


Ex(p)r © 1 Ox are components of the plas- 
i yp) ¥(p) E1lC Stain 


Dr From Hooke's Law 


Ey(e) = 1/E(6y = V6,) 
Oxy (e) = 1/6 (Txy) 
where E 1S Young's Modulus 
Gis the Shear Modulus 


Y/Y is Poisson's Ratio 


c. Stresses must Satisfy equilibrium conditions: 
O6y/ax + dtxy/dY 
Otky/ ax + O6y/oY 


qd. Total strain must satisfy compatibility: 


Letting E' = Ze a 9 %¢ 2 
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and R' | + Zc ee 
ayes © A Syees © Sse ies 
Then E' + R' = 0 by compatibility. 
If R' is non-zero (uneven plastic strain distribution), 
then elastic (residual) stresses must exist to satisfy 
compatibility. Implications of the above are: 
(1) If plastic strains are linear with respect to dis- 
placement, no residual stressesS occur. 
(2) Residual stresses cannot be determined by measuring 
the stress change during external loading or unloading. 
(3) Residual stresses 6,y,6y, Jxy can be calculated if 
Ex(e)s Ey(e)s Xxy(e) are determined. 
The residual stress iS measured by meaSuring the change 
in elastic strain after cutting a specimen. Residual 
stresses can then be calculated uSing the following 


equations: 


6x = -E/(1-42) (Ey + Ey) 


6 


y = Sy Ae + wEy) 
where the above strains are experimentally measured. 

By measuring the Strains, at the same location, on both top 
and bottom of the specimen, one can determine the approximate 
contribution of planar and bending stresses to the total 
residual stress. The average of the strains through the 
plate will give the planar stresses. The difference in 
Strains through the plate will give the contribution from 


bending. 
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3.3 Experimental Procedure 

Twenty strain gages were attached to the Sine curve 
fabricated as one of the experiments in reference [2]. A 
photograph and parameters used to make it 1S Shown in 
Appendix A. Ten strain gagesS were placed on both sides at 
the same x-y coordinate. The plate is a 12" x 12" x 1/2" 
mild steel plate. Since the primary deflection is perpendi- 
cular to laser beam path (parallel to the y axis), the focus 
of the experiment was determining the residual stress in the 
x direction. Two gages per side were placed parallel to the 
beam path to determine if a "poisson" effect caused strain in 
the y direction. Figure 3.3-1 shows the plate geometry and 


gage placement. 


y 10 GAGES PER SIDE 





Figure 3.3-l: Placement of Strain Gages on the Sine Shape 
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All gages were BLH Electronics SR-4 type FAE-25-35-S6EL. 
baechenaemed gage faccor of 2.02 + 1%, a resistance of 350 + .5 
ohms. K had a value of +0.26 and the gage length was 6.35 mm 
(1/4 in) long. 

A Vishay Instruments Model P-350A Digital Strain 
Indicator with two SB-l Switch and Balance units was used to 
measure all Strains. This instrument allows the measurement 
of ten strain gageS per Switch unit. One-half, full, or one- 
Quarter balanced bridge circuits can be selected. For this 
experiment, the one-quarter bridge configuration was used. 
In this mode temperature effects can be automatically compen- 
Sated by connecting an extra lead from one of the strain gage 
connections to an internai dummy resistor (350 ohm) within 
the strain indicator. A 5000 micro-Strain (pin/in) internal 
calibrator and a gage factor input circuit allows the compen- 
Sation of the measured resistances for the wire resistance. 

Each strain gage is balanced to "Zero" strain indication 
with the Switch and Balance unit prior to cutting the speci- 
men. After cutting the specimen, the main indicator control 
is varied to rebalance the bridge circuit and the strain in 
microinches-per-inch is read directly off the Strain Indica- 
tor unit. Table 3.3-1 shows the major specifications of the 


maicator Unit [7]. 
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Range feo ,U00 maicrormch/inch polar- 
ity reversible. 


Sensitivity Variable. Null meter deflects 
from zero to full scale with 
40-4000 pin/in at gage factor 


e 


Readability lpin/in. 


Accuracy mule or redding or 5 pin/7in 
whichever 1S greater for R=120 
ohms, GF=2. +0.3% of reading 
or 5 win/in, whichever is 
greater for R=120 ohms, GF=1.5 


orr4 . 

Gage Factor Continuously variable from 0.1 
POmel 0. Ol. 

Bridge Excitation 1.5 VRMS at 1000 HZ Square wave. 


Table 3.3-l: Vishay Instrument Strain Indicator and Balance 
and Switch Unit Specifications 
3.4 Experimental Results 
Table 3.4-1 shows the results of strain measurements of 
the sine shape. Longitudinal spacing is relative to the y 
axis. Unless otherwise specified all gages are transversly 
placed relative to the y axis. MicroStrain differences and 


averages are between the top and bottom of the specimen. 
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DIST. FROM MICROSTRAIN MICROSTRAIN MICROSTRAIN MICROSTRAIN 


CENTERLINE TOP BOTTOM DIFFERRENCE AVERAGE 
-4 tla +9 () +30 +106 
-3(LONG) -~182 -394 -~212  -288 
-3(TRANS) -582 -458 -1040 -620 
-2 -560 -620 -60 -~590 
-3/8 -245 -~468 -~223 -357 
+3/8 = 25 -~440 -~222 -~329 
+2 = eo.6 -~288 -~152 -~212 
+3 (LONG) -462 =o -~156 -~540 
+3 (TRANS) -~1498 +1359 -2857 -695 
+4 -20 -180 -160 -100 


Table 3.4-1: Experimental MeaSurements of Micro-Strain for 
the Sine Shape Formed by Laser Bending 


At point x = -3 the ratio of strain in the longitudinal 
to transverse direction is: -182/-582 = 0.31. At point x = 
wo, the ratio is: -462/-1498 = 0.308. Both values are 
approximately poisson's ratio for steel (V4 = 0.3) showing 
that bending of this plate in the direction parallel to the 
laser pasSs iS primarily a "poisson" effect. This concept 
Will be extremely important inthe analyses of the deflec- 
tions of the complex shapes and the subsequent development of 
the algorithms for them. Furthermore, on the basis of these 
results, the values of the elastic strain in the y direction 
at all points on the plate will be assumed to be: Ey(e) = 


Ex(e) during the residual stress calculations. 
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The residual stresses due to bending and planar effects 


now can be calculated using the following equations: 


6, = -E/(1-¥) (Ey + VEY) 
—- aa 2 & 

Sy = -E/(1-¥*) (Ey + VE) 

Ex =¥%€, (based on experimental evidence) 
For steel: E = 2.9 x 10? pPpSl and Y= 0.3. The equations 
Become; 

by = -6.37 xX 10+(€,) Wieles Eas the measured 

6, = -1.16 x 107(€y) strain in pin/in. 


Table 3.4-2 shows the residual stresses calculated from the 
experimental data given in Table 3.4-1l1 uSing the above 


equations. Figure 3.4-1 shows the results graphically. 


Bret. FROM MICROSTRAIN MICROSTRAIN X-STRESS X-STRESS 
CENTERLINE DIFP ERENCE AVERAGE BENDING PLANAR 
(KSI) (KST) 
~4 5.0 106 ~1.91 ~6.76 
ao = AG Ia one 18.4 
-2 = 610 -590 Sao S726 
a / a8 — Shee) 14.2 22.8 
3/6 mae ~329 14.2 210 
I 2 ore = NMG oro See 
+ 3 eG ~540 o.29 4 34.4 
+4 -~160 -~100 G2 Ga) 


Table 3.4-2: Bending and Planar Residual Stresses Determined 
for the Sine Shape Formed by Laser Bending 
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The bending stress at x = #2 to the edge is 
approximately constant since the edge is clamped on that 
edge. The bending stress reaches a local minimum and the 
planar stress reaches a maximum at approximately the position 
of maximum magnitude of deflection where the rate of change 
of the deflection iS a Minimum. Therefore, for this shape, 
the bending residual stress iS a minimum where the slope of 
the deflection (dz/dx) is a minimum. Likewise, the planar 
residual stress iS a maximum where the magnitude of the 
deflection (|z|) is a maximum. Except for the region of 
external clamping of the plate, the bending stress is 
consistently less than the planar stress. The mean ratio 
between bending stress and planar stress for the positive 
Momecet1On halft-cycle (xX < U0) is 0.43 and the mean ratio 
between bending stress and total residual stress is 
approximately 0.48. The fact that the bending residual 
stress represents about one-half of the total residual stress 
could indicate that significant changes in the magnitude of 
the deflection might occur if stress relieving is performed 
on the plate. Future experiments could confirm this 
hypothesis. 

3.5 Summary and Conclusions 
The results of the measurement of residual stresses in 


the sine shape appear to indicate the following: 
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fee sralhi simaguecea in the plate ina direction 
parallel to the laser pass is approximately poisson's ratio 
times the strain in the direction perpendicular to the 
heating path. 

b. The bending residual stress 1S a minimum where the 
deflections are a maximum since the magnitude of the slope of 
the deflection is a minimum at that point. The planar resi- 
dual stress iS a maximum where the magnitude of the deflec- 
tions are a maximum. 

Ce The bending residual stress constitutes, on the 
average, of about 50% of the total residual stress but varies 
on the plate from about 10% to about 803%. 

qd. Although clamping the plate appears to make uniform 
the bending stress distribution, clamping does not appear to 
affect significantly the planar stress more than about one 
beam diameter (~1.5 inch) away from the restraint. 

e. The planar stress is approximately symetrical about 
the centerline of the plate less the effect from artificially 
restraining one side, 

fF. No local variance was visible in the planar stress 
due to multiple laser passes. This indicates that for over- 
lapping or close proximity laser line passes, the stress is 
distributed about the plate in accordance with the geometry 
of the shape. 

g. The highest y direction residual stress was observed 


to be about the yield strength of the metal (765 ksi). 
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CHAPTER FOUR 


PHASE II LASER LINE HEATING EXPERIMENTS 


4.1 Introduction 


This set of experiments was part of the second phase in 


the sequence whose objective is the development of a closed- 


loop controlled process for laser bending of steel] plates. 


The specific goals of these experiments, conducted at Naval 


Research Laboratories, were: 


a. 


Verification of and testing the reproducability of 
previously obtained data. 

Increasing the data base by investigating the behav- 
HOMmmOr 35/8 andu5/6sineh thick»=mild steel plate. 
Testing a simple linear algorithm for predicting the 
deflections for simple curves as a function of plate 
thickness and heat delivered to the metals surface. 
Quantitatively Investigating the effect of multiple 
laser passes at different locations. 

Investigating the effect of surface preparation on 
the deflections through multiple pass-single loca- 


tion experiments. 


As discussed in Chapter 2, reference [1] discusses the 


Get lectionmof stee® plates as awfunction of the heat 


delivered to the metals surface. The power delivered to the 


Sok 





metal's surface is represented by the following parameters: 
P//T or P/t/V 

where P is the laser power delivered to the metal, t is the 

plate thickness, and V is the bedplate/steel speed beneath 

the laser. 

It 1S anticipated that even if the process of laser 
bending steel plates to produce complex shapes 1s controlled 
by a closed-loop system, initial bending to within twenty to 
thirty per cent of the final shape will be governed by the 
data base and/or developed algorithms. Final shaping and 
dewarping would be monitored by the optical laser/video 
System. A major goal of these experiments was to teSt a 
Simple, linear algorithm to predict the one dimensional 
deflection of plates with thicknesses not previously within 
the data base. To examine the reproducability of this 
process, the deflections of 1/2 inch thick plate were 
measured and compared to that obtained during Phase I experi- 
ments. 

Since the completed system/process may require multiple 
passes at the same or different locations on the plate, both 
multiple pass-single location (MPSL) and single pass-multiple 
location (SPML) experiments were conducted. It was expected 
that fundamental question would be answered concerning any 
interaction of passes from one location to another and inter- 


actions of multiple passes at the same location. 
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4.2 Experimental Procedure 

In order to minimize the number of Variables and to 
Simulate a probable industrial environment, generated laser 
power was held constant at 7 kW with a beam diameter of about 
1.5 inches. Beam spot geometry was the "top hat" discussed 
in reference [2]. Variation of the heat linear density 
delivered to the metal was accomplished by varying the plate 
Speed past the Stationary laser beam. No complex shapes were 
created in this series of experiments. Beam path locations 
were predetermined on the basis of: 

a. Maximizing beam pass Separation when independent 
measurementsS were desired and controlling the separation 
distance during SPML studies. 

IF Obtaining no effects from the plate's edge parallel 
Co the travel direction, 

(oe Allowable clearances on the bedplate mechanism with 

the steel plate, 

All experiments were performed with 24" x 36" K-TEN20CF 
steel plates. Each waS spray painted black. No clamping or 
forced cooldown was done. Deflections were measured with 
three dial gages of which two were placed on one side of the 
beam pass, Deflections were taken Poorodically during 
cooling until no change was observed (about 20 minutes per 
pass). Figure 4.2-1 shows the arrangement for a typical 


experiment. 
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DiAL_ GAGES 





Figure 4.2-l: Typical Arrangement for Phase II Experiments 
Conducted at N.R.L. (September 1985) 


4.3 Experimental Results 

The three dial gage readings were plotted and, assuming 
the minimum deflection to be at the centerline of the pass, 
tme angular deflection on both sides of the pass was 


determined from: 


1,2 = |e - 2 ey 
X12 “Kp 
Of = Oj + 09 
where 212 are the measured deflections at coordinate 


X1 98 Xp is the coordinate of the centerline of 
’ 


the pass, and O¢ is the final total angular deflection. 
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For an ideal plate, the minimum deflection at the centerline 
of the beam pass, 21,2min, should be Zero. In relality, 
deflections parallel to the beam path are superimposed on the 
plate resulting in a uniform (except at the edges) raising of 
the plate. This is probably due to the longitudinal residual 
stresses discussed in Chapter Three. Since only x direction 
deflections are desired, this overall plate translation in 
the z direction was removed by subtracting the minimum value 
which was found to be consistently greater than zero. The 
algorithms for complex shapes developed in the next chapter 
will include this effect. 

The linear algorithm assumes that S6¢ is a linear fun- 
ction of P/W between any 2 adjacent, previously determined 
data points. Furthermore, it assumes S8¢ is a linear function 
of the thickness for a given P/W between any two adjacent 
thicknesses. Deflection angles for 3/8" thick plate was 
peocted using values of 6, for 1/4" and 3/8" plates from 
figure 2.2-1 or table 2.2-1]1. 1/2" thick plate was subjected 
to values of P/VV not pecyvTousciy vwithane the data base al- 
though they are shown graphically in figure 2.2-1. Values of 
Of were predicted on the basis of linearity between 
bracketing values of O¢ within the data base (table 2.2-1). 
This was done to test the linearity of the curves within 
neighborhoods of the preexisting data and to test the repro- 


Giicamim tty Of “Ene process. Finally, 5/8" thick plate was 
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subjected to three different P/\V and the deflection angles 
measured. Table 4.3-1 gives the results of these experiments 
and the predicted values for 3/8" and 1/2" plates. Figures 
4.3-1 and 4.3~2 are graphs of these meaSurementsS Superimposed 


upon the results from Phase I experiments. 


THICKNESS P/N P/t\V 6 8 
ACTUAL PRED. 
nie 3.41 0.28 0.50 
1.49 3.98 0.86 1.19 
3/8" 
1.75 4.67 1.23 eo 
2.21 5.90 1.66 125 
1.49 2.98 0.47 0.42 
2" 1.75 3.50 0.73 0.70 
2a 4.43 1.04 OW 
1.49 2.39 0.38 
5/8" 1.75 2.80 0.82 
7 3.54 0.98 


Table 4.3-l1: Test Results Obtained in Parametric Studies 
Conducted September 1985 
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On the basis of the experimental resultsS, one can 
observe: 

a. The close agreement between the previous and current 
results for the 1/2" plate indicating the reproducability of 
the process. 

De For 1/2" plating, a linear algorithm predicting 
deflections between data points previously obtained is 
reasonably accurate. 

c. 1/4" plating should not be used to linearly predict 
the behavior of plates with a different thickness. 

a. The measured deflection angles for the 5/8" plating 
were consistently higher than expected with or without the 
use of a linear algorithm. PO: S ureason, —1t 15S 
recommended that additional experiments be performed on 5/8" 
plating prior to the integration of these results into the 
data base. 

Experiments were performed to determine the effect of 
multiple passes on the same plate. P//V was held constant at 
1.278 and the plate thickness was 3/8". The results are 
given in table 4.3-2. The cumulative (total) deflection 
angle as a function of the number of passes is shown in 
figure 4.3-3. The results show that for as little separation 
as one beam diameter (71.5 inches) between passes the 
deflection 1S approximately constant. This implies that 
under the conditions of this experiment, superposition of 


angular deflection is a reasonable assumption in SPMNL mode. 
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PASS DIST. TO DIST. TO e, 
NEAR. PASS NEAR. END (degrees) 
(inches) (inches) 
A : 10.5 24 
B Does 5 oe 
C as 10.8 Be 
D eae Gas 38 
E 3.0 13.5 ee 
F 2.88 7.88 31 
G 1.50 es 33 


Average Angular Deflection Per Pass 18S 0.3 Degrees 


Table 4.3-2: The Effect of Pass Separation Distance on the 
Angular Deflection of 3/8" plate with P/\V = 1.28 
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Total Angular Deflection (degrees) 


ia 


1.9 
1.8 
ie 
1.6 
1.5 
1.4 
1.3 
IBZ 
1.1 


0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 


Figure 4.5-3: Number of Passes vs. 


Angular Deflection (t=3/8", P//V=1.28) 


5 | . 


Number of Passes 
— §= .30N - .08 
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The last series of test in this sequence indicated that, 
without surface preparation (painting) after each pass in the 
MPSL mode, no Significant deflection occurs after the initial 
pass. It should be emphasized that in these experiments the 
plate was allowed to cool to ambience after each pass. Since 
the absorption of ~10 micron radiation by a metal increases 
rapidly with increasing plate temperature, multiple passes 
at a Single location might indeed cauSe significant 
deflection if the material is not allowed to cool 
significantly prior to the next pass. Further investigations 
of the effects of surface preparation are recommended. This 
should include, at a minimum, the use of alternative 
preparations such as white titanium oxide which is frequently 
used in laser scattering experiments and should include the 
MPSL processing without cooling after each pass. 

4.4 Summary 

The results of the series of Phase II experiments condu- 
cted at Naval Research LaboratorieS in WaShington, D.C. con- 
tribute to the understanding of the behavior of metals when 
subjected to laser line heating. They also provided valuable 
data in the development of a closed loop control system. 
Specifically: 

a. Results for 1/2" plating were verified showing the 
reproducability of the process. 

b. Multiply located passes on the same plate cause 


deflections which are approximately linear with respect to 
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the number of passes. This was confirmed using K-TEN20CF 
with beam path separation to as low as one beam diameter, 

c. If no surface preparation is performed and the plate 
is allowed to cool to room temperature after each pass inthe 
MPSL mode, no meaSurable deflection occurs after the first 
pass. 

a.. mice Orme nhm 1S SUEFIiclient to predict 
deflections for a given thickness at power densities not 
Within the data base. A linear algorithm using 1/4" data for 
predicting behavior of plates of other thicknesses is 
inaccurate. 

Sufficient predictability and control has been shown to give 
encouragement that development of a closed loop laser bending 


system 1S practicable for simple shapes. 
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CHAPTER FIVE 
AN ALGORITHMIC APPROACH TO PREDICTING DEFLECTIONS 


IN A COMPLEX SHAPE 


5.1 Introduction 

Very elaborate techniques have been developed or applied 
to the prediction of the behavior of metals undergoing line 
heating. In spite of this, to date, there exists no algo- 
rithm based on either engineering physics or a data base 
Which can synthesize a complex, three dimensional shape. 
Recalling the discussion in Chapter Two about the overall 
goals of present research in laser bending, algorithms will 
be essential for automating a shape synthesizing process 
uSing line heat. This is also probably true for automating a 
metal straightening (dewarping) process. This is equally 
true for laser and flame line heating. 

The "purest" algorithms are those which are based on a 
thorough understanding of the physics of the technique 
including the mechanisms of behavior of the object in ques- 
Earn . This approach is usually possible only for simple 
systems where there are few primary variables. Even then, 
the use of very powerful computers or experimental equipment 
is often required. Any algorithmic based or assisted auto- 
mated technique requiring the use of a main frame computer 
Will probably not be cost effective for a manufacturer such 


aS a shipbuilder. 
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The other major approach is to base the algorithm on 
empiracal data. This approach is often preferable to the 
"Dhysics" approach since it allows the development of a 
process with potentially many primary variables and an incom- 
plete knowledge of the science of the phenomenon. A main 
disadvantage of this approach is that it may be more expen- 
Sive to develop the algorithm because of the number of exper- 
iments required to create a sufficiently large data base. In 
the long run, for producing an automated process for shaping 
steel by thermomechanical means, the later approach is best 
Since it will reduce the complexity and cost of the hardware 
in the complete system. In fact, it is not unreasonable to 
Strive for a system whose control system is maintained and 
driven by a microcomputer. Therefore, the algorithms should 
be aS simple and compact as possible. 

The main objective of this development is not to produce 
the most elaborate, comprehensive, or sophisticated algorithm 
for laser bending steel plates into complex shapes. Rather, 
it is to prove that simple linearized algorithms can be 
constructed which will produce with reasonable accuracy those 
Shapes. Furthermore, these algorithms must be based on pre- 
existing information on simple, Single pass deflections al- 
ready within the data base. Although the data base for 
Simple angular deflections is substantial, the number of 
preexisting complex shapes is not. Therefore, it is fully 


expected that, aS more experiments are performed on complex 
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Shapes and aS more analysis is performed on the mechanisms of 
bending in response to line heating, the algorithms will be 
modified to become more general and more accurate. This 
paper develops one algorithm for each of the shapes 
considered. These are the parabaloid (dish), the parabolic 
hyperboloid (saddle), the cone, and the antisymmetric double 
paraboloid ("sine"). 
5.2 General Methodology 

The algorithms created in this chapter are empirically 
based on only one sample of that particular shape. Although 
a strenuous attempt was made to include known physical pheno- 
mena, on occasion, underivable factors had to be included 
Which may or may not remain the same for the same shapes 
constructed under different conditions. The same methodology 
waS applied to the different shapes whenever possible. 
Linearity was assumed whenever excesSSive loss of accuracy did 
not result. Symmetry was determined from geometric proper- 
ties of the shapes. A photograph of each of the analysed 
shapes is included in Appendix A. These are extracted from 
reference [2]. Also included with the photograph is a sketch 
of the pattern of laser passes and relevant parameters used 
to construct the shape. The deflections at specific points 
on the shape were measured with a dial gage. All data was 
compensated for the thickness of the plate and was normalized 
such that the deflection at the center of the plate (x,y = 0) 


waS numerically zero. The coordinate system is cartesian 
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With the center of the plate the center of the coordinate 
system. The width of the plate in the "x" direction is 
notated W, and the width (height) of the plate in the "y" 
direction is notated Wy: Unless otherwise noted, the 
"top" of the plate is that side such that the laser path is 
in the +y to -y direction. For a perfectly continuous pure 
geometric shape, the entire shape can be generated by func- 
tions for the x = 0 and y = 0 contours. Figure 5.2-1 shows 
how the saddle shape can be generated with only the contours 
of x = 0 and y = 0 analytically determined. Throughout this 
Chapter "z(x,y)" will represent the deflection at point (x,y) 
and z(0,y) and z(x,0) will represent the components of the 


deflection from the x = 0 and y = 0 contours, respectively. 


Zi 





Z(OY) =gfY) AX =X) + glY¥)* CONSTANT 


Figure 5.2-1: 3 Dimensional Shape Generated from Contour Lines 
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The total deflection at a point (x,y) is: 
Zx, yy) = £(x) | yao + gly) |y=o + Constant 
Since the curve is normalized such that 2(0,0) = 0 the 


above constant iS zero and: 


z(x,y) = £(X) | yig + Gly) lyan oF 
ex, ¥) = ex, Oy + 200; Y). 
Geometric symmetry simplifies algorithm generation. For 


example, geometric symmetry predicts for the perfect 


bowl and saddle: Z(x,y) = z(-x,y) = z(x,-y) = z(-x,-y); 
For the sine shape: Z(x,y) = -z(-x,y): 
For the cone shape: Z(x,y) = Z(-x,y) only. 


The above relationships are valid only if the center of the 
coordinate system is at the center of the shape as is assumed 
for all the analysis herein. 

The actual complex shapes were prepared in such a way as 
to try to maintain the symmetry. This was done by the choice 
of beam pattern and the equality of laser irradiation. 
However, when the shapes were analyzed it was found that the 
Symmetry was not always preserved. This could be due to 
unevenesS in Surface preparation or plate thickness, 
irregularities in laSer power or bed plate speed, or unequal 
edge effects. These deviations of the real plates geometry 
from the ideal geometric shape can be categorized as: 

a. Warping: Lack of fairness or Smoothness for a 
particular z-contour line. 


b. Asymmetry along an entire contour line or surface. 
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Figure 5.3.l1(a)-l, which will be discussed ina later sec- 
tion, shows examples of both types of irregularity. 
Geometry and the pattern of laser passes should cause us to 
predict that for a given y = k line the deflection at +x is 
equal to the deflection at -x. However, the deflections in 
the y>0, x>0 quadrant appear consistently greater that the 
values in the y>0, x<0 quadrant. Warpage is also evident for 
the y = 5 line and from x = -3 to x = -4 on the y = 1 line. 
For all the shapes the edge effects are somewhat erratic. 
The dish was irradiated in such a way to minimize edge 
effects. Nonetheless, it still appears in several of the 
contour lines. For some of the other shapes not fabricated 
to minimize end effects the edge effectsS may or may not be 
considerable. 

Due to these irregularities, the following general 
methodology was used to analyze the shapes: 

a. The raw (actual) deflections were normalized such 
that z2(0,0) = Q. 

b. USing a fairing process with splines commonly used 
in naval architecture to smoothe hull shapes, the contour 
lines were "dewarped". This process was performed such that 
the deflection was the same value at a point (n,k) approached 
from both a y = k and a x =n contour. 

c. Deflections were averaged in such a way as to obtain 
geometric symetry predicted from the ideal curve. 


qd. The dewarped, symmetrized shape, called the model 


49 





Shape, (2,), was then used to generate the algorithm (Z,). 

Although for each of the Shapes the algorithmically 
determined deflections (Z,) are compared with the measured 
values, the goal wasS to construct the modeled (dewarped and 


symmetrized) deflections (2 to within about 10 % accuracy. 


m? 
This procedure can be justified by the fact that plate 
warpage and lack of complete symmetry always cannot be 
controlled or predicted uSing a real process. Secondly, 
these algorithms are designed to initialize the automated 
process. That is, these algorithms should be capable of 
obtaining "in the ball park" deflections. It is assumed that 
the closed loop feedback system will be capable of fine 
tuning the shape which would include dewarping and 
symmetrizing the shape as required by the application of the 
Shape. It is estimated at this time that approximately 70 to 
85 %$ of the desired shape's deflection will be obtained by 
use of the algorithm and/or the data base. The remaining 
deflection will be controlled by the feedback system. 

The process of constructing the model shape consists of 
using the single pass @6- (aS a function of P/\V) from 
reference [2], the number and pattern of laser passes, 
generalized effects from the edges, near-pass coupling 
effects and, as required, correlating factors to generate the 
y = 0 and x = 0 contours. These contours are then used to 
predict the deflection at any given point within the plate in 


the manner described above. 
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The slopes of the x = 0 and y = 0 model contours were 
mot constant. In most of the shapes the contours were 
divided into three regimes with their own governing 
equations. These regimes are the flattened area near the 
maximum or minimum, the edge and corresponding nearby region, 
and that portion between the other two regimes. In several 
cases the last two regimes had the same slope and were then 
considered to have two regimes. TO improve the accuracy near 
the origin where flattening of the shape was prominent, a 
Single point deflection is predicted half-way between the 
Origin and the borderline between regime one and two. 

Some method had to be developed which would correlate 
the angular deflection within the data base for a single 
laser pass, the number of passes on that portion of the plate 
and a factor accounting for beam path coupling effects and 
edge effects which were distributed over all the plate. 
Defining "k" to be the distributed angular deflection: 


] where 


kK = n*0,-°[Wy (1+2) °B.D. JEW, - (14+7)°B.D. 
2°Wy Wy 


nis the number of laser passes on that line or plane of 
symmetry; 
Se is the angular deflection (Figure or Table 2.2-1); 


W., and W,, are the widths of the plate; 


x ye 
B.D. 1S the beam diameter; 
Y is poisson's ratio. 


For steel with the laser beam diameter of 1.5, (1+t7/)B.D. is 


approximately equal to 2. The equation becomes: 


om 





K =en°*O,°(W. - 2) (W.. - 2) 
2W,W 


iy 


The factor of two above comes from the definition of @,¢, 
The distributed deflection slope 1s given by: 

Az/\Kq = Tan(k) where q is the generalized 
coordinate. 

For example, a1/2" thick plate with width of 24" and 
height of Poe is symmetrically irradiated with 20 passes on 
one side to produce a dish. The beam power is 6.3 kW 
including optical losses and bed plate/steel speed is 12 ipm. 
Then p/\V = 2.02 and from table 2.2-1, re = eo 0. From the 
above equation: 


k (20) (0.9) (22) (16) /(2) (24) (18) = 7.339 


Tan(k) = Tan(7.339) = 0.1286 = Az/{q 
The distributed deflection angle is the basis of the algo- 
rithms of all the shapes. 2Z,(x,0) and z,(0,y) are assumed to 
be functions of the Tan(k) parameter. It should be stated at 
this time that the Tan(k) parameter is more accurate than 
Superposition (n°@-/2) for three dimensional curvature with 
many passes. However, for simple one dimensional passes such 
as those described in chapter four and references [1] and 


(2], Superposition appears to be more accurate. 
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5.3 Shape Analysis and Algorithm Generation 
5.3.1 Dish Shape 
a. Derivation of the Model and Algorithm 

The dish shape was produced by irradiating a plate on 
one side longitudinally, transversely, and diagonally to the 
y = 0 coordinate line. This pattern is shown in Appendix A 
with relevant plate and heating parameters. Table 5.3.1-1 
gives the measured and normalized (raw) deflections and 
figures 5.3.1-1 to 4 are graphs of the measured data. 

Using the fairing and symmetriZzing procedure described 
in the previous section, model contours are produced. Table 
5.3.1-2 and figure 5.3.1-5 show the model deflections. 

Recall that the distributed deflection slope is given by: 
Az/Aq = tan(k) = tan[n@¢ (Wy-2) (W,-2)/2W,Wy]} 


Berethis dish: n = 30; Wy=W =10"; O¢ =0.7 which gives 


By 
k = 6.72° and tan(k) = tan 6.72° = 0.1178 

Since the geometry of the dish which was used to generate the 

model has x/y symmetry, the y = 0 and x = 0 contours are the 

Same curve. Using this fact and the value of tan(k) above, 


the following equations can be used to construct the model 


dish. 
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Figure 5.35.1—5: Model Deflections 


After Fine] Symmetrizing 


MODEL LINEAR DEFLECTIONS 





oY=0 +1Yi=1 o M=2 RI Qs x M=4 J iyis5 


X y=+/-5 Y=+/-4 Y=+/-3 Y=+/-2 Y=+/-1 Y=0 
Eve Pees DE Eoin. DE eee DEP iy. Dine. DE 2 ees 
1/5 se 22 ay el 601 72 Ud ~433 0 
+/-4 ea) e602): ~492 oo 22.0 ou 
+/-3 647 465 342 243 pales eas 
ct = 2 oo ioe ~- 248 ~-148 ~ 086 Os 
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Table 5.3.1-2: Model Deflections After Final Symmetrizing 
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Za(x,0) = (#//2)tan(k) |x| x = 0.1W, 

Za(O,y) = W/2)tan(k)ly| y = 0.1Wy 

Za(x,0) = YWtan(k) |x| Orel, as PX 0. ew, 
Z,(0O,y) = Ytan(k) |y| O.1W, < ly| < 0.2Wy 
Za(x,0) = tan(k) [|x| - .2W, (1-1 ] .2W, < |x| < .4W, 
ZalO,y) = tan(k)Ely| - .2W,(17/)} -2Wy < |y| < .4Wy 


Za(x,0) = tan(k)[.75|x| + .2Wy-.1lW,] .4Wy, < [x] < .5W 


Xx 


Za(O,y) = tan(k)(.75|y| + .2Wy-.1lWy] .4W 


a 
A 
a 

[A 


e Wy 


The multiplicative factor of .75 in the slope of the last 
series of equations accounts for edge effects. Figure 5.3.l1- 
6 shows the algorithmic and model deflections for |q| = 0. 
To find the deflection at any point on the plate one now only 
has to add Z2,(x,0) and z,(0,y) for the particular coordinate 
(x,y). 
Or Comparison of the Algorithm, Model, and Measured 
Deflections 

Table 5.3.1-3 gives the model and algorithmic 
deflections and the deviations in per cent between them. The 


average deviation between zZ, and z, is about 8.5% which is 


a 
within the goal of 10%. Outside one beam diameter radius of 
the center where flattening is greatest, the average 
deviation is about 7.5%. Figures 5.3.1-7 and 5.3.1-8 show 


the predicted and model curves for the dish shape. Figures 


5.3.1-9 and 5.3.1-10 show the predicted and measured curves. 
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Pagunecans.o..-1. amd 5.3.1-12 show a “3-D" perspective of the 
dish shape. Figure 5.3.1-11 shows the measured deflections 
While figure 5.3.1-12 show a 3 times magnification of the 


measured deflections to better show the Shape and warpage. 
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Table 5.3.1-3: Model and Algorithm Deflections 
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X Y Z mdl Z alg DEVIATION 


or /-} 7 a % 
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Table 5.3.1-3: Model and Algorithm Deflections (continued) 
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2.3.2 Cone Shape 
a. Derivation of the Model and Algorithm 

The pattern of laser passes and relevant parameters used 
to produce the cone shape are given in Appendix A. Table 
5.3.2-l1 gives the actual, normalized data. Figures 5.3.2-1 
to 3 are graphs of the raw data. In this case, imposed y- 
direction strains are present due to both the poisson effect 
and due to the non-orthogonality of the beam passes relative 
to the coordinate system. Therefore, longitudinal bending 
has a contribution approximately proportional to the sine of 
the angle between the y axis and the beam pass. 

Analysis of the data seemed to indicate a significant 
coupling effect between beam passes. The x-direction 
deflections are significantly greater (2.5-5.75) than 
expected by either superposition (n@¢/2) or by the tan(k) 
parameter. This significant deviation from expected behavior 
may be an inertial effect. As the angle of the beam pass 
with the y axis increases, the closer the beam is to the 
corner of the plate. The much less inertia of the plate on 
the outboard side of the beam pass should imply that the 
majority of the angular deflection from the beam pass would 
be manifested on that outboard side. As the angle increases, 
the larger this effect would be. Another possible 
contributor to this "amplification" effect is the density of 
heating at the focus of the beam pattern. Geometric symmetry 


dictates that, for the model, 2(x,y) = z2(-x,y). 
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No symmetry exists across the y = 0 line. Table 5.3.2-2 and 
figure 5.3.2-4 and 5.3.2-5 give the symmetrized and dewarped 
model deflections. AS one would predict from the geometry of 
a cone, the changes in deflection slopes as a function of x 
are significantly greater for the region of Y > 0 than for 
the region of Y < 0. For this reason one algorithm for each 
region was developed. Since the plate was cut as shown in 
Appendix A, the effective area term in the tan(k) parameter 
was formulated on the basis of the area enclosed by the 
perimeter of the shape minus 2 inches divided by the area of 
the plate. Using these values gives: 
tan(k) = tan[(7)(.7)(282)/(2) (356.25)] = .0337 

It was discussed above that an "amplification" factor was 
evident in the values of the deflection slopes of the model. 
One set of empiracal equations that satisfy the observed 
increases in deflection slope are: 


f = ] x < W tan(aj,) + BD 


Y 


a= Wy Sin(a,)/k W tan(a,)+BD < x < W 


y, - y 
f= Wysin(a3)/2°k’*Y Wytan(a3) + BD < x < W,/2 


tan (a3) + BD 


y: 


where f is the amplification factor; aj yedio7 23) 4a5e the 
angles of the beam pass with the y axis; kis the distributed 
deflection angle; BD is the beam diameter; and vis poisson's 


ratio, For the conditions under which this cone was made: 
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Table 5.3.2-2: Model Deflections for the Cone Shape 
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Directly 


equations, the algorithm's equations become: 
Z,(0,y) 


Z,(0,y) 


I. y = 0 


Z_(x,0) 
Z_(x,0) 


Za (x,0) 


ire y < 0 


Za(x,0) 
Z_(x,0) 


Z_(x,0) 


mit. y > 0 
Z_(x,0) 
Z_(x,0) 


Z,(x,0)=tan(k)(5.75(1- 


cE nS) Soe SS Wenz 


-(1-7Y)tan(k) [ly] - 3] 


2.5tan(k)[x - 0.1W,) 


Cant) I> yx — 4W/ 3] 


2,oban(k) 1x ——0. lwy) 


eeW Ry (SoS I Ee = Aas Ie) 


inputting the above values for f into the 


-Wy/6 < ¥ < Wy/6 
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NN 
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NN 
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OE XG 16 


Wy/6 < x < W,/3 


tan(k) [(1-Y) (5.75) x-.758W,}] 9 W,/3<x<W,/2 


0 << xk W,/6 


~—_ 


Wy /6<x<Wy/3 


Y_) - 1.1375W, +11.5K%) 


Z(x,0) = tan(k){5.75x-1.1375W, + 11.57} 


geaex, 0) = Z.(-x,0) 


for x 


Figure 5.3.2-6 and 5.3.2-7 are graphs of the model and 


algorithm y = 


0 and x 


0 contours. 





Ol 8 Vv O 


LO 
Z'O 
SO 
%°0 
fo S°0 
9° 
£°0 
270 


INOZUOD O=A YD BUOD JOL BSUOI}DIIJOqG eD 


WUIYP4OB|\y puo jepoW :9—T' EG s4Nnb!l4 


NOLWWOS 14450 YVSaNN 


8 2 





O=|plui--x os O=Bj0— es 
(NID A es 


— ~ 


ANOJUGD O=X }O Bucy 4oj SUdIVOaIJEqG 


GL‘O— 
vL‘O— 
eL°O— 
ZL°O— 
LL°O— 
L“O— 
60°0— 
g0°0— 
L0°0— 
90°0— 
so’O— 
vo°'O— 
£0°0— 
ZO°O— 
Lo‘O— 
O 


WUWOB|\yY pub jepoW :7—Te'G s4anbl4 


NOOSA WSC YVSNIT 


8 3 





De Comparison of the Algorithm, Model, and Measured 
Deflections 

Table 5.3.2-3 shows the predicted and model deflections 
and the deviation between the two. The average deviation is 
5.2%. Figures 5.3.2-8 and 5.3.2-9 compare the predicted and 
model curves. Figures 5.3.2-10 and 5.3.2-l11l1 compare the 
predicted and measured values. Figure 5.3.2-12 shows a "3-D" 


perspective of the measured values. 
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Table 5.3.2-3: Model and Algorithm Deflections 
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Table 5.3.2-3: Model and Algorithm Deflections (continued) 
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5.3.3 Sine Shape 
a. Derivation of the Model and Algorithm 
Appendix A shows the parameters and the pattern used to 


form the "Sine" Shape. The coordinate system used in the 


development of the algorithm is 1809 from that shown in the 
Photograph in the appendix. Table 5.3.3-1l1 gives the 
norlmalized, raw data. Figures 5.3.3-1 to 5.3.3-4 show the 
raw data. As can be seen from the graphs, the edge along the 
x = +6 line was clamped. Geometric symmetry dictates that 
Meee yY) = -z(X,yY) and 2(x,y) = 2z(x,-y). Deflections in a 
direction parallel to the beam passes were found. Their 
Origin has been discussed previously. The clamping 
complicates the modeling and algorithm formation Since it 
causes violation of the symmetry. Since all the shapes have 
and will assume unrestrained plates only the half-plate which 
is unclamped will be used for symmetrizing. 

The first attempt to model the plate focussed to fitting 
a Sine curve to the shapes. This proved to be too 
mraccurate. Instead, the shape wasS modeled as an 
antisymmetric double parabola. Taolkema.S.s—Zzeana Ligure 
3.5.3-5 shows the modeled data. The inflection point is 
approximately at x,y =0. Significant coupling exists between 
the two sides of the plate. The result is that the portion 
of the plate containing the line connecting the two sides of 
the plate (each side irradiated on opposite sides) has an 


unexpected rate of change of curvature. ThisS area actS aS a 
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clamp. Therefore, for this reason, the tan(k) parameter must 
be modified by poisson's ratio from about x = -BD to about 
Xx = +BD. The tan(k) parameter remains applicable to most of 
the interior of the plate. 

The following equations form the heart of the algorithm 
for the sine shape: 


tan(k) = tan[(4) (0.7) (100)/(2) (144)) = .0172 


Z,(O,y) = -.25tan(k) ly] ly| < Wy/6 
Za(O,y) = -tan(k)[]y] - 15/Wy] Wy/6 < ly| < Wy/2 
Z,(0,y) = 2,(0,-y) < 0 


For the "negative" deflection half-cycle: 


Z,(x,0) = -(1//) tan(k) x O< x < Wy/6 
Za(x,0) = -tan(k) [x - Wy(%-1)/6°*27] Wy/6<x<.94W,/4 
Za(x,0) = tan(k)[x - .846W,] 94W,/4 < x <Uwy 
Z,(x,0) = tan(k)[x - W,/2] UYwy < x < Wy,/2 
ae T.1e7 ; ; p 


For all regions: 
Ze OZ CX, 0) 

b. Comparison of the Algorithm, Model, and Measured 
Deflections 

Table 5.3.3-3 gives the results of the algorithm and the 
model with their deviations. Note that this algorithm is 
somewhat inaccurate for computing the values at the edges. 
This is in part due to the clamping effect discussed earlier. 
The average deviation more than one beam diameter from the 


edges is 9.2%. Figures 5.3.3-6 and 5.3.3-7 show the model 
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and algorithm deflections for x,y = 0. Figure 5.3.3-8 shows 
the model and algorithm deflections for x,y not equal to 
zero. The +/- y contour deflections for a given x are the 
same value for both the algorithm and model. Figures 5.3.3-9 
and 5.3.3-10 show a comparison of the actual and algorithm 
deflections. Since the x = +6 edge is clamped which distorts 
the actual deflections, only the x < 0 deflections are 
compared. Figures 5.3.3-ll and 5.3.3-12 are 3-D plots of the 


Sine shape. 
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Figure 


3-D Perspective of the Sine Shape with 


Actual Deflections 
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Figure 5.3.3-12: 3-D Perspective of the Sine Shape with 


Deflections Magnified by 3 Times 
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5.3.4 Saddle Shape 
a. Derivation of the Model and Algorithm 

Appendix A shows the parameters and pattern in which the 
plate was fabricated. Table 5.3.4-1 and figures 5.3.4-1 to 
5.3.4-4 show the actual normalized deflections. For this 
shape, symmetry dictates that z(x,y) = z(-x,y) and z(x,y) = 
z(x,-y). USing the above rules of symmetry and the general 
methodology utilyzed for the other shapes, the model deflec- 
tion data was determined and is given in table 5.3.4-2 and 
shown in figure 5.3.4-5. 

Ideally, since the power delivered to each side of the 
plate is the same, the magnitudes of the deflections should 
be the same. This is not the case. The magnitudes of the 
deflections of the X = 0 contour are much more than predicted 
from the tan(k) parameter or superposition. It is believed 
that the poisson effect which has been discussed is 
responsible. The "top" plate irradiation imposes longitudi- 
nal deflections which enhance the deflections during the 
"bottom" plate irradiation. Therefore, the side which has 
been irradiated second will have more deflection than 
predicted by superposition or tan(k) for single side 
irradiation. This increase should be predictable from pois- 
son's ratio and, in fact, the plate can be reasonably modeled 
over much of the area by assuming a multiplicative (ole) 
factor with tan(k). It appears from the data that the second 


Side passes do not reduce the first side (y = 0 contour) 
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deflections which are reasonably modeled by the tan(k) 
parameter. Perhaps, the reason is that the process has 
insufficient power to overcome the high longitudinal residual 
stresses imposed during the "top" plate irradiation. On the 
other hand no such stresses oppose the longitudinal bending 
Guring the first series of passes. 


The algorithm equations become: 


tan(k) = tan[(16)(0.7)(.64)/2}] =0.0626 
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EpeeompPparison of the Algorithm, Model, and Measured 
Deflections 

Table 5.3.4-3 shows the predicted and model deflections. 
Average deviation is 10.6%. The average deviation on the 
plate greater than a distance of about 2/3 BD from the edges 
is 9%. Figures 5.3.4-6 to 5.3.4-8 compare the results of the 
algorithm and the model. Figures 5.3.4-9 and 5.3.4-10 
compare the results of the algorithm with the actual 
deflections. Figures 5.3.4-ll and 5.3.4-12 are 3-D 


perspectives of the saddle shape. 
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Table 5.3.4-3: Model and Algorithm Deflections for the Saddle 
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Table 5.3.4-3: Model & Algorithm Defl. for the Saddle (cont) 
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Figure 5.3.4-ll: 3-D Perspective of the Saddle Shape with 


Actual Deflections 
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3-D Perspective of the Saddle Shape with 
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CHAPTER SIX 


CONCLUSION AND RECOMMENDATIONS FOR FUTURE RESEARCH 


The main purposes of the research and analyses inthis 
paper have been proven. Namely: 

a. The technique of irradiating a plate with a laser 
under controlled conditions is a reproducable and predictable 
phenomenon. Complex shapes can be mathematically constructed 
with reasonable accuracy uSing simple, linearized equations 
and well-known aspects of geometric symmetry. 

On Generality between algorithms has been established 
through the use of the TAN(k) parameter, poisson's ratio, 
plate dimensions, and the same general methodology. 

c. Residual stress analysis appears to be consistent 
with the general algorithmic development with respect to 
poisson's effect. 

qd. All methods used to establish the construction algo- 
rithms are applicable to flame line heating as well as laser 
line heating with appropriate modifications for differences 
in heating parameters. This should provide sufficient flexi- 
bility for further development of a general line heating 
automated process using either lasers or flames. 

It was stated previously that the algorithms developed 
in this paper are only a "first step". Admittedly, they lose 


accuracy near the edges and sometimes near the points of 


eG 





maximum and minimum curvature. Future research should, in 
part, focus toward "fine-tuning" them. Futhermore, the even- 
tual goal of the algorithmic development should be to genera- 
lize the shapes into as few algorithms as is practicable. 
Possibly all conical and spherical shapes can be generalized 
into a single algorithm. Due to the lack of experimental 
data, it was not possible to generate general algorithms 
which would cover a variety of pure or combined shapes. For 
example, no algorithm was generated which would predict with 
correlatable accuracy a combined shape such aS a sine curve 
Superimposed on a conical shape. Also, Since no experimental 
data exists on shapes which were unsymetrically heated in 
terms of the P/ V for the various passes, the applicability 
of these algorithms to this condition is untried and may not 
be accurate. However, futher experimentation including that 
which is on-going at the time of this writing should clarify 
and better define the scope of applicability of these 
routines. 

Most research, up to this time, has concentrated in 
producing a process to laser bend metals without much atten- 
tion on the mechanisms and consequences of the process. The 
practicability of the process has been demonstrated. At some 
point the physics of laser line heating should be thoroughly 
investigated. This is of more than academic importance. For 
example, current studies indicate that laser beam spot geome- 


try has a significant effect on the magnitude of deflections. 


ed 





Beebo turns Out tO be tue, P/ V, in itself, may no longer 
be the best parameter to uSe to correlate deflection angles. 
More research into the questions of power, thermal, and 
stress distributions are recommended. 

With the exception of the conical shape, all the complex 
Shapes formed up to now have been fabricated with beam passes 
either parallel or perpendicular to the edges of the plate. 
This may not be the most efficient pattern. Other patterns 
to produce the same shape may reduce the warpage and 
asymmetry found in the shapes already constructed. A dish 
might be formed by a circular or spiralling beam pattern and 
it conceivably could be more efficient and/or produce less 
warpage. Therefore, research is recommended in studying 
Various patterns to produce a shape and in correlating the 
total heat required to produce the shape as compared to the 
patterns used in the experiments described in this paper. 

Perhaps the most important of the recommendations is the 
investigation of using the laser bending with closed loop 
System to symetrize and dewarp a plate made with the laser. 
It has been shown that almost any shape can be fabricated. 
What has not been demonstrated up to now is if the same laser 
system can economically produce those shapes to the Same 
tolerances required in many industries including ship 
bila Ldang. The use of the laSer interferometer/video system 


Will be invaluable in accomplishing this. 


2 





REFERENCES 


1. Masubuchi, K., "Phase I Report on Laser Forming of Steel 
Pores £Or Ship Construction", submitted to Todd Shipyards 


Corporation, Los Angelese Division, February 1985. 


2. M°Carthy, R. W., "Thermomechanical Forming of Steel Plates 
Using Laser Line Heat", S.M. Thesis, Department of Ocean 


Engineering, M.I.T., Cambridge, MA, June 1985. 


3. Deacon, D. L., Material Degradation in Heavy Plates Caused 
by Bending with a Laser", S.M. Thesis, Department of Ocean 


Engineering, M.I.T., Cambridge, MA., June 1984. 


4. Haidemenopoulos, G. N., "Heat Flow and Material 
Degradation During Laser Metal Forming", S.M. Thesis, 
Department of Ocean Engineering, M.I.T., Cambridge, MA., June 


IG 5. 


>. Moshaiov, A., "Thermo-Elastic-Plastic Plate Bending by a 
Boundary Element Method with Initial Plastic Moments", Ph.D 
Dissertation, Department of Naval Architecture and Marine 


Engineering, University of Michigan, Ann Arbor, Mich., 1985. 


6. Masubuchi, K., Analysis of Welded Structures, Oxford: 


Pergamon Press, 1980. 


7. Vishay Instruments,"Model P-350A Digital Strain Indicator 


Instruction Manual", November 1975. 


133 





Appendix A 


Photographs: Complex Curvature Study 
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-HEAT INPUT 

BEAM DIAMETER 
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3/3 inch 
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6.3 KW 

13 ipm 

31.5 KJ/in 

1.5 in (38.1 mm) 
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CONE SEGMENT 


PLATE SIZE 18in X 34in 
PLATE THICKNESS 1/3 inch 
MATERIAL A-36 Mild Stee! 
POWER INPUT 6.0 KW 

PLATE TRAVEL SPEED 8 ipm 

HEAT INPUT 37.6 KJ/in 

BEAM DIAMETER 18/8 in (43.1 mm) 


NUMBER OF PASSES 7 (topside only) 
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SOLID LINE: TOPSIDE PASS 
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SINE SHAPE 


PLATE SIZE 13in X 13in 
PLATE THICKNESS 1/3 inch 
MATERIAL A-26 Mild Stee! 
POWER INPUT 6.3 KW 

PLATE TRAVEL SPEED 13 Ipm 

HEAT INPUT 31.5 KJ/in 
BEAM DIAMETER 1.5 in (38.1 mm) 
NUMBER OF PASSES 8 (4 each aide) 








SOLID LINE: TOPSIDE PASS 
DASHED LINE: SIGE OPPOSITE 


We s i: 





SINE SHAPE 
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PLATE SIZE 

PLATE THICKNESS 
MATERIAL 

POWER INPUT 

PLATE TRAVEL SPEED 
HEAT INPUT 

BEAM DIAMETER 


NUMBER OF PASSES 


SADDLE SHAPE 


10in X 10in 

1/3 ineh 

A-38 Mild Stee! 
6.3 KW 

13 ipm 

31.5 KJ/in 

1.6 in (38.1 mm) 


32 (16 long. top, 16 tran: 
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SADOLE SHAPE 


SOLID LINE: TOPSIDE PASS 
DASHED — OPPOSITE 


vq {" 


147 





po) renee re aes 


Wi 
ME Eas 
ue es 


a rit 


yg 


a oe , 
ie ee - 3 


Oe ein ‘3 or Pe be ole bf a Naf te Py H DY 


a 
ane. Fi 
ae ry en : La 

yy Viet 


tHe an 
A re Pet 
ee 


esate Few Hsee 


eye wes 
etn fo ce La ri 


iP AE 
" a%. a eee he 
a: y : i eT Le 
Ce Mi: am a 
RP ie 1,% 
| re a Sy’ 
@ f a 


, te 
hae ay 


ue 


4 : ‘ 4 eat i . os 7 sf i ji - Tae ry 
ve 


ret Rn ‘ta . “™%]. ’ ¢ , , nt Se J s peat rad 


" oa Sa A 4 is 
. bee» a) Dy D rs ad A 
il We aera ney ; Ny ae tah a a 


aa is PIR al si > wer age te 4 : ; ‘ : ft. A 


oe A : a G = F ae 
are | be 
0 ; m oe 


ca 
CO ee aA 
a ly 





148 














t 





| cy i aie = Or, 
7% aha. 


Pig hn Atel tig. hd Mh 6 Ae <a ” 
pie Sl LR i AAA Oe res ae ; one 
y, Q t weit ' SS ARE hie Th . Ae | oA SH 
aah af Bel dats Paine eles Time Y 
SAREE GUAM Xi pd es, 


yA? . 4. : 

daw be Be es, 

bs 2 eal) . . } ff ae Pr A 5 4 ' - = 

te oS wie ; al Catan oat: Secon diay Ty eae Me ey rae ar 
j : i iz , o : ; ul —- ? 4 ac a 








| or 




















> 
5 

















Pi they EES eg wey 
ech Ue Madey ix ’ i 
pelae aes ees : 4 ARS : " 


on 
a myst 8 + 
i 








" 
, ye 





ek wh iaay 
; Sees eat a 

. ie ae dae, : 4 15 ate 7 ' 
Mr At Gi Lie 





iy ds .d * res mt 41 ' k 
OS NTA Hea 
tel 





at Ae sine 


er th ee 
Tah ae ay 





lay 





APPENDIX B 


Application of the Algorithms/Predictions 


for Future Experiments 
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Ongoing research at the United Technologies Research 
Semucr VU LRG) imemiartford, Connecticut by M.I.T. staff and 
Students is based on increaSing the number of samples of 
complex shapes and, using the algorithms within this paper, 
to predict the deflections under given power parameters. En- 
closed within Appendix B are the predictions for the 
deflections (SsSymmetrized and dewarped) of the shapes. 
Furthermore, the research will determine if "fine tuning" 
With a laser is practicable. That is, can the formed plates 
be irradiated again to achieve the final desired deflections? 
All plates within this appendix are assumed to be 
18"x24"xl/2" and the pattern of laser passes are the same as 
for those plates used to develop the algorithms. P/V and 
the number of laser passes are varied to determine the 
accuracy of the algorithm. 

The beam spot of the CO» laser at UTRC is annular vice 
the "top hat" spot geometry of the laser at NRL. Beam 
diameter was varied in accordance with optical power theory 
to give approximately the same power density as that 
Gelivered by the 6 KW laser at NRL. The accuracy of this 
approach will be determined only after the analysis of the 


samples. 


oe 





I. DISH 


A. Parameters 


-P/W 1.708 
~O¢ Get = 
-~Number of Passes 10 diagonal 


10 diagonal 
10 longitudinal 
10 transverse 


B. Equations 


tan(k) = tan[(40)(.417) (22)916)/(2) (24) (18) ] 
zZ(x,0) = 0.0179[x| | ieee 5 " 

z(0,y) = 0.0179[y| yal S62" 

Zi, Oy = 0 0retetea| | 2. Sime. |e 5" 
Z(0,y) = 0.0358|/y| Dean Y | ano... 
z(x,0) = .1192|x| - .4005 Son X | Ga on 
z(0,y) = .1192|y| - .3004 3.6" < ly] < 7.2" 
z(x,0) = .0894/x] - .1144 9.6" < |x] < 12" 
z(0,y) = .0894|y| - .0858 7.2" < ly| < 9" 


C. Results 
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II. DISH 


A. Parameters 


=e 2S 

4 Dee 

-Number of Passes Same as I. 

B. Equations 

lO) ees 1S 
z(x,0) = 0.0533]/x| Pxcleoue2 5" 
Zi 077 = 0.0533) y) eva 
z(x,0) = 0.1066|x| Deo || <5)! 
z(0,y) = 0.1066|y| De orl yl en SG 
z(x,0) = .3553|x] - 1.1938 niles 9.6" 
z(y,0) = .3553|y| - .8954 On) | Vale eee 
Z( Xo) = 26.6.5 x | Said Ora | Xoo 
ZAOG Yeas 2605 | Y lo Zoe Bis 2 |'¥ |e" 


C. Results 
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III. CONE 


A. Parameters 





=P/ WV 1.708 

-OF 0.417° 

-Number of Passes Same aS Appendix A 

B. Equations 

z(0,y) = 0 ie SS 
Z(0O,y) = -0.0140[y - 3] ese 
zZ(0,y) = 0.0140[y + 3] y < -3" 
Y = 0 
z(x,0) = .0201|x| oes 
z(x,0) = .0501[|x] - 2.4] 4" < [x{ < g" ** 
Z(x,0) = .1153|/x| - .6416 hy ema ox |e 2 a 


*This equation applies for all y's (+/- or 0) 


*"This equation applies for Y < 0 








Y < 0 
z(x,0) = .0807[x| - .3649 Speer ieee 
Y > 0 
Z(x,0) = .0737|[x| - .2145 Cael 
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C. 


Results 


eZ 


12 


12 


- 180 


396 


~742 


060 
~ 228 
~456 


~670 


~060 


noe 


04! 


7oOs 


-.042 


.018 
- 186 
~407 


-628 


12 


dg2 


1D 


~.042 


~018 


oie 


Lg 


220 


-.084 


-.024 


~143 


~ 364 


7205 


-.084 


-~.024 


096 
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IV. 


CONE 


A. 


B. 


Parameter 


Equations 
ZO, y) 


zZ(0,y) 


2 


o 


ZX 0) 
ZA x0) 


Ze) 


Ss 
-P//V 
-~O¢ 


-Number of Passes 


= 0 

= -0.0404[|y| - 1.8485] 
=e Oo cn xa 

= 0.1444[|x| - 2.4] 
Sees 2 Se ee) Sloe 


Same as Appendix 


| Y | 


ly | 


g" 


[A 


3° 


< 


S 


*This equation applies for all y's (+/- or 


**This equation applies for Y < 0 


=). 23 Sights 0 Sea)) 


.2123[|x{ - .6181) 
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C. Results 
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12 


eZ 


eZ 


LZ 


ey 3 
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2.14 


hela 


Uz 


eZ 


12 


eZ 
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e121 


~-052 


~ 399 


~920 


-~.243 


-.070 


~413 
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Sine 


A. Parameters 


-P/WV 225 
-O¢ ig Dae 
-Number of Passes 8 Passes Each Side 


B. Equations 


tan(k) = tan[{(8)(1.2)(.815)/2] = 0.0684 


Z(x,0) = ~0.2279x x < 4" 

Z(x,0) = -.0684x - .6383 LSS SESS Ce 
z(x,0) = .0684x - 1.3885 Sa BCX << 7 oe 
ZU Oy = .193 exe 2.5108 eee tee Ceres 
Z(x,0) = -z(-x,0) 

ZO) = = 05129 /[y | YAS" 
z(0,y) = -.0684|/y| + .0570 Aan 71 eno) 


iRexs 





C. Results 


X Me Z X Y Z 
0 0 0 0 6 eS OS 
3 0 -.684 3 6 =-1,04 
6 0 -.978 6 6 ec 
9 0 era 2 6 Seles 
eZ 0 = LS) S V2 6 -.546 
aS 0 ~-684 3 6 ~331 
0 0 UES =0 6 OZ 
=e 0 es = 6 ~ 420 
ame 0 oS eZ 6 = OO 
0 3 -.154 0 9 = 5 Ses 
3 Ss 7630 3 9 —1,24 
6 3 = eee 6 2 -1.54 
9 3 eee 2 9 ore 
eZ 3 -.34/7 eZ 9 -./51 
= 3 2530 ao 9 0 
-6 3 ~824 =a0 9 ~420 
-9 3 »619 -9 9 Zo 
ely Z 3 »039 ee 2 =o Or 


oe 





VI. SINE 


A. Parameters 


-P/ NV 


-OF 


-Number of Passes 


Bae tiaGiens 


ZG, oe= —-. 0791x 

Z(x,0) = -.0237x - .7383 
Z(x,0) = .0237x - .4818 
ZAXGO) =. COTOXe— 2.8713 
Z(x,0) = -2(-x,y) 

z(0,y) = -.0178|y| 
ZAG) = —.0237 
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eae 


0.417° 


Same as VI. 


a. eX 


1A 
Uy 


SrA. Xx 


[A 


oe Sex 


ly] < 4" 


an < fy| < 


{A 


Woe. 
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Cc. Results 


Xx vi Z X b Z 
0 0 0 0 6 —serleZ S 
3 0 ee 3 3 6 Oe 
6 0 =e S 6 6 -.462 
9 0 =~. 268 9 6 ie Sel 
eZ 0 ao EZ 6 Se eho, 
3 0 es = 3 6 ~114 
= 6 0 soa9 = 6 - 216 
-9 0 - 268 -9 6 ~145 
ale 0 Oy =k 7 6 nora 
0 3 es 0 2g -.194 
3 s = 7290 3 9 -.431 
6 Ss =r 2 6 9 =o aes 
9 3 erie 9 g -.462 
eZ 3 =. 120 MZ 3 = Ore 
= 3 5 -184 =3 9 ~043 
0 3 seb = 9 ~145 
ao 3 72> -9 9 .074 
bz 3 ~014 aoe 9 rl 2.) 
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VII. SADDLE 


A. Parameters 


B. Equations 


tan(k) 


ZX 0) 


ZUOr yy} 


ZC x7 0) 


Z (0, Vv) 


ZX 5.0!) 


Zr0,y) 


ZO Ry } 


where: 


-P//V 


-~OF¢ 


-Number of Passes 


tan[(16)(.417)(.815)/2] = 


-008 


Ons 
.0158|x| 


.0260|/y| 


.0475|x| 


-0617|y| 


.0867ly| 


ZX Oe 0; 


ea 


0.41 


8 


ae 


16 Top; 16 Bottom 


0475 

ere 

Ne es ee 

Dee << el ee 15 

2" < ly] < 3.6" 
1520 50 (xj) < alt 
.1286 3.6" < ly| < 7.2" 
~3084 eee aya eo 


(O52) Mae 
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C. Results 


X Y Z X x Z 
0 0 0 0 6 -.242 
3 0 -048 5 6 ~wl 94 
6 0 ers 6 6 =O 
9 0 a2) 9 6 5031S 
12 0 ~420 IL 7 6 iS 
0 3 seme 0 o -.472 
3 B ree 3 9 -.424 
6 3 Ore 6 9 ore 
9 3 ok oy 9 9 = ome 
We 3 ~ 342 12 9 arent Eee 
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VIII. 


a. 


Bi 


SADDLE 


Parameter 


Equations 


tan (kK ) 


a(x, 0) 


z2(0,y) 


z(x,0) 


Z (Oey) 


z(x,0) 


2VO cy) 


z(0,Y) 


where: 


Ss 


-P/ WV 


-OF¢ 


-Number of Passes 


Cen( (16) (1.2) ¢.815) /2) 


#0229 


-0376 


.0458/x| 


07 Sistly| 


isla) - « 


-LPEGiy| - . 


.2509|y| - . 


PX a0) > O; 


Same as VII. 


= bsive 
;e) 2.5" 
Fy) =oe8 


Zoe <a < 


a < iy| < 3.6" 


fea 


4397 Ba <j ah << 72" 

7 oy Sao" =< tive < 722" 
e623 fee a 
Z(0O,y) < 0. 
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VIII. SADDLE 


A. Parameters 


-O¢ Mee 
-Number of Passes Same as VII. 


B. Equations 


Comme mtamitlo) (1.2) (.815)/2) = .1374 

Z (x0) = .0229 | ms 5" 

20g ae. 037 6 vie= 2 

ZC xO 7. 014 58 | x | oe <x) <5" 
ZO = 20753) y | ely | <13e6." 
z2(x,0) = .1374|[x| - .4397 eens eee 2 
z(0,y) = .1786|y| - .3721 3.6" < |y| < 7.2" 
z(0,y) = .2509|y| - .8923 pee lay < Olt 


Where: 2(x,0) > 0; 2(0,y) < OQ. 
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